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Abstract 



We study electromagnetic waves scattering by a 2D photonic crystal made of 

a stack of diffraction gratings. In case where there are only two propagative 

modes in the crystal, we derive an explicit expression for the superior (resp. 

inferior) envelope of the reflected (resp. transmitted) energy. 
Photonic crystals are periodically modulated devices having photonic band gaps The 
theoretical study of these structures may be performed by use of Bloch waves theory, multi- 
scattering theory or they can be modelled as stacks of gratings which are the 
elementary layers of the device (cf. fig. 1). We use this last approach in this note. 
Two situations can be roughly distinguished in the scattering behaviour of periodic struc- 
tures: when the wavelength is of the same order of the period, the resonant domain, there 
are infinitely many evanescent modes and finitely many propagative modes in the structure 
(the number of which being given by Bloch waves theory). As the wavelength increases, 
the number of propagative modes diminishes and is finally reduced to 2. In the very large 
wavelength regime, the structure behaves as if it were homogeneous. In case of stacks of 
gratings, that leads to replace one grating by a homogeneous slab (fig.l) with some effective 
permittivity depending on the polarization [^,0. However, there is an intermediate zone of 
wavelengths for which there are only two propagative modes inside the gratings, but the 
wavelength is not large enough to use a homogenization scheme (for instance A is only two 
or three times the period). This work deals with this last situation. The crystal is made as 
a stack of N gratings (fig 1) with vacuum above and below. 

We consider first an elementary layer. As there are only two propagative waves, there exists 

a 2 X 2 matrix Ta,^ such that T f ^_ j = f j (see fig. 2 for the notations), the elements 

of T are denoted (tj^). This matrix can be computed rigorously from a standard theory of 
diffraction by gratings (for instance the Coupled Waves Method, see § for example). The 
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main property of matrix T is that det (T) = 1 [|T0| so that if /i is an eigenvalue then fi ^ is 



the other eigenvalue. We have shown in another article that in a layer characterized by a 
2x2 transfer matrix, the reflected and transmitted coefficients for layers can be obtained 
in close form : 

u^^ - 1 

rNiX,0) = f ^ 



- fg-' 

where /i is an eigenvalue of T associated with an eigenvector {ui,U2), an eigenvector associ- 
ated to is denoted by (^1,^2) and 

f ^ = -J5 — [ — . 9 A, e) = — ■ 

ifJoVi + V2 iPoUl + U2 (2) 

and = — cos{d) 
A 

Remarks: 

1) It is worth noticing that functions / and g do not depend on a particular choice of an 
eigenvector. 

2) The expressions are valid whatever the angle of incidence and polarization, it is 
only the particular expression of matrix T that changes. 

The expressions given in (|1|) show that it suffices to compute the transfer matrix of 
one layer to compute the transmitted and reflected fields. However, the reflected 
and transmitted energies are highly oscillating functions which can be better represented by 
their envelope (superior and inferior limits when N tends to infinity). We have shown in |Tl[ 



that the superior envelope R^o of the refiected energy, and conversely the inferior envelope 
of the transmitted energy Too are given by : 



T 



A-tr (T? 



(ti2/5o - t2i/?o"') (3) 

Rcia 1 -Ty 
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We now turn to some numerical examples. We use as an elementary layer a lamellar dielectric 
grating with an air layer below (see fig. 1). The parameters are the following: h/X = 1, Sa = 
4 and £6 = 1 with a filling factor of 1/2. It is illuminated by an s-polarized monochromatic 
plane wave under normal incidence. We have plotted in fig. 3 the refiected energy versus 
X/d, obtained for a stack of 20 layers. We see that above Xc/d = 1.515 the true reflection 
ratio is well approximated by (|1|): this corresponds to the critical wavelength above which 
there are only two modes in the structure (cf. fig. 4 for the band structure). Above Xc/d we 
see an almost perfect agreement between the refiection computed from (|1]) and the refiection 
obtained from a rigorous code (fig. 3). Finally, we have plotted the envelope of the refiected 
energy obtained from (^, we see a perfect agreement above wavelength Xc/d. 
We have developed an approximate theory for the study of photonic band gap materials in 
case where there are two propagative modes in the structure in both s- and p-polarization. 
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This method allows explicit computations of the field once the 2x2 transfer matrix of one 
layer is known. This last matrix can be obtained numerically, then equations (|l]) and @ 
can be used to obtain close-form expressions of the diffracted field. This can be useful for 
studying the so-called "ultra-refractive" effects ||T^-[14 . 
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Figure captions 



Figure 1: Geometry of the photonic structure : the elementary cell is made of a grating, 
with period d filled with materials whose relative permittivities are and Eb- There is a 
dielectric layer below. 

Figure 2: Definition of the transfer matrix of a layer. 

Figure 3: Reflection vs X/d. Solid lines : reflection obtained from a rigorous theory. 
Dots : reflection obtained from formulas (|ip. 

Figure 4: Band structure of the photonic device: when A > Ac there are only two 
propagative waves in the crystal. 

Figure 5: (a) Envelope of the reflected energy (bold line) obtained from (|^) and reflected 
energy (solid line). 

(b) Zoom of figure 5 (a) around the bifurcation value Ac/ d. 
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Figure 1 




Figure 2 



